Correlated Random Walks and the Joint Survival 

Probability 



Mark B. Wise a and Vineer Bhansali 6 

(a) California Institute of Technology, Pasadena CA 91125 
wise@theory . caltech . edu 

(b) PIMCO, 840 Newport Center Drive, Suite 300 
Newport Beach, CA 92660 
bhansaliOpimco . com 



Abstract 

First passage models, where corporate assets undergo correlated random walks and 
a company defaults if its assets fall below a threshold provide an attractive framework 
for modeling the default process. Typical one year default correlations are small, i.e., 
of order a few percent, but nonetheless including correlations is very important, for 
managing portfolio credit risk and pricing some credit derivatives (e.g. first to default 
baskets). In first passage models the exact dependence of the joint survival probability 
of more than two firms on their asset correlations is not known. We derive an expression 
for the dependence of the joint survival probability of n firms on their asset correlations 
using first order perturbation theory in the correlations. It includes all terms that 
are linear in the correlations but neglects effects of quadratic and higher order. For 
constant time independent correlations we compare the first passage model expression 
for the joint survival probability with what a multivariate normal Copula function gives. 
As a practical application of our results we calculate the dependence of the five year 
joint survival probability for five basic industrials on their asset correlations. 



1 Introduction 



In the management of credit risk and the pricing of corporate bonds and credit derivatives 
a fundamental role is played by corporate default and survival probabilities. Structural first 
passage models [1-7] are one of the standard approaches to estimating these probabilities. 
In such models the firm asset value evolves with time in a random walk and default occurs if 
the value of the firm's assets fall below a default threshold. Merton [1] showed that corporate 
stock and debt are options on the firm asset value. Equity holders hold a call option on 
the firm asset value, while bond holders are short a put on the firm asset value. In first 
passage models a simple analytic expression can be derived for the risk neutral probability 
of a company surviving without defaulting during a time interval t, P(t) 1 . Risk neutral 
probabilities are used for pricing corporate bonds and credit derivatives. A no arbitrage 
argument implies they are calculated by taking the drift in the firm value to be given by 
the risk free rate of return rather than its "real world" value. 

For portfolios of bonds the joint survival probabilities are needed for a proper assessment 
portfolio risk. In addition the pricing of some credit derivatives (e.g., first to default baskets) 
depends on the risk neutral joint survival probabilities. These can be calculated using a 
first passage model with correlated random walks or by postulating a particular Copula 
function. A multivariate normal Copula function is often used. In first passage models 
the exact dependence of the joint survival probability of more than two firms on their 
asset correlations is not known. We derive a formula for the joint survival probability 
that includes the effects of asset correlations treating the correlations to linear order in 
perturbation theory. For constant time independent correlations we compare this expression 
for joint survival probabilities with what a multivariate normal Copula function gives. At 
linear order in perturbation theory in the asset correlations the joint survival probability of 
n firms can be expressed in terms of the firms default correlations. 

The perturbation theory results derived in this paper for the joint survival probability 
are valid provided the asset correlations and the number of firms are not too large. We 
argue that for five firms perturbation theory provides a useful approximation to the joint 
survival probability for realistic values of the asset correlations. When the number of firms 
is very large our results provide a limit (i.e., the small correlation limit) where numerical 
simulations of correlated random walks, with first passage boundary conditions, can be 
compared with analytic results. 

As a practical application of our results we calculate the dependence of the five year risk 
neutral joint survival probability of five basic industrials on their asset correlations. The 
companies we consider are: Alcoa Inc., Dow Chemical Company, E.I. du Pont de Nemours 
and Company, International Paper Company, and Weyerhauser Company. 

2 Joint Survival Probabilities 

The asset value V\ of a company % is the sum of the value of its stock Si(Vi, t) and its debt 
Bi(Vi,t). Assuming dividend payments to stock holders of company i, qiVidt in a time 
interval dt, we take the firm asset values of the n-firms (i.e., i = 1, . . . ,ra) to evolve with 

lr The risk neutral probability of the company defaulting during the time t is then 1 — P(t). 
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time according to 

n 

dVi = mVidt + Vi ^2 Vijdwj - qiVidt, (1) 
j'=i 

where Wi are standard Brownian motions, i.e., 

E[dwidwj] = dtSij (2) 

and the Kronecker delta is defined by 5ij = 1 if i = j and Sij = if i ^ j. In Eq.(l) ^ is 
the drift of the z'th company value if there are no dividend payments to its stock holders 2 . 
Here we are assuming that dividend payments are financed out of operating profits or cash 
reserves and not out of the issuance of new debt or stock. We also assume that the dividend 
yields qi are constant independent of time. Time depended dividend payments were studied 
in [8]. The assumptions made in this section are consistent with the Modigliani- Miller 
theorem [9,10]. 

Eqs. (1) and (2) imply that the asset (return) covariance matrix pij is 

n 

Pij = a ik a jk- (3) 
k=l 

Hence the asset volatilities are 



<Ti = 



a X! a ik a ik, (4) 
\ k=l 



and the asset correlation matrix is 



1 n 

iij = — ^2 a ik a ik- (5) 



UiU i fc =i 



We allow the default barrier to have an exponential time dependence and assume that 
company i defaults at a time t if its firm value falls below, 

V id = di e^, (6) 

where di and Aj are constants independent of time. This is the framework of Black and Cox 
[2]. One usually takes the exponential dependence to be given by the risk free rate r, i.e., 
Aj = r. This form for could arise from a covenant that firm i has with its debt holders 
which forces the firm into bankruptcy to ensure that the debt holders do not lose more than 
a fixed fraction of the present value of their principal. 

The default probabilities follow from the expressions for the evolution of the firm value 
and the default barrier given in Eqs. (1) and (6). The constant real world asset drift fii in 
Eq. (1) does not affect £j(V$,t) and Bi(Vi,t). Risk neutral default and survival probabilities 
that are used for pricing corporate bonds and credit default swaps are obtained by setting 
the asset drift ^ equal to the risk free rate r. 



2 To simplify the discussion in this section, we take the debt to be zero coupon bonds that mature at a 
date that is beyond our investment time horizon. 
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It is convenient to change variables from the firm value Vi to Zi defined by, 

Zi{Vi,t) = QogiVi/Vio) - V). (7) 
Default occurs if Zi falls below the time independent default threshold 

z id = log(d i /V i0 ). (8) 

In Eq. (7) Vio is the initial asset value for firm i. Ito's lemma implies that the time evolution 
of Zi is given by, 

n 

dzi = rjidt + ^2 Vijdwj, (9) 
i=i 

where 

m = m,- ^ - h - qi (io) 

Vi can be expressed in terms of Z{ and t using Eq. (7). 

Let p(z±, ... ,z„,, t)dzi • • • dz n be the joint probability that the n companies survive to 
time t (i.e., none of them default before time t) with each Zi in the interval dz{. More 
explicitly, 

p(zi, ... ,z n , t)dz 1 ■ ■ ■ dz n = P[tdefauit > t for all n firms|,Zj(i) G dz^. (11) 

The probability that all n companies survive to time t (i.e., the joint survival probability) 
is given by 

Pi,...,n(t) = dz 1 --- dz n p(zi,...,z n ,t). (12) 

Since Zi evolves via Eq (9) the probability density p(z\, . . . , z n , t) obeys the diffusion 
differential equation (see, for example, [11]), 

l + E^-^E ^8^j - » W) = 0. (13) 

where the sums go over all n values of i and j. The probability density p(zi, . . . , z n , t) obeys 
the initial condition 3 (here Zi > z^) 

p(z 1 ,...,z n ,0) =S(z 1 )---S(z n ), (14) 

and the first passage boundary conditions, 

p{z ld ,z 2 ,...z n ,t) = ... =p{zi,...,z n - 1 ,z nd ,t) = 
p(oo,z 2 , ■ ■ ■ z n ,t) = ... =p(z 1 ,...,z n - 1 ,oo,t) = 0. (15) 

It is convenient to reexpress this in terms of the volatilities and the correlation matrix, 
d d 1 ^ 2 d 2 1 ^ _ d 2 \ , 

at + - 2 ?^ a? - 2 ^'^ = °- (16) 



3 S(x) denotes the Dirac delta function. It is denned by, 5(x) = for x =fc and J S(x)f(x) = /(0) for any 
smooth function f(x). 
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For n = 2 the solution to this differential equation with the appropriate boundary conditions 
for the first passage problem has been expressed in terms of an infinite series [12]. Here we 
consider the case where the off diagonal elements of the correlation matrix are small enough 
to treat them as a perturbation. Expanding the joint survival probability density in powers 
of the off diagonal correlation matrix elements we write 



p(z 1 ,...,z n ,t) =p(°\z 1 ,...,z n ,t) +p {1) (z 1 ,...,z n ,t) + ... 



(17) 



where p(°\zi, . . . , z n ,t) is the solution with the asset correlations set to zero, p^\zi, . . . ,z n ,t) 
contains all the terms linear in the off diagonal elements of the correlation matrix, etc. A 
similar expansion is made for the joint survival probability itself 



A,..,»W = ^ ( l„(t) + ^.,„(*) + ---- 



(1) 



(18) 



We consider the case where the volatilities and correlations are constants independent of 
time. Then 

P { °\zi, ...,Z n ,t) =p 1 (z 1 ,t) ■ ■■p n (z n ,t), (19) 

where the probability density pj(zj,t) is defined so that pj(zj,t)dzj is the probability that 
the company j survives to time t (i.e., does not default before time t) with Zj in the interval 
dzj. It is given by 



p(zj,t) 



2<KG 2 t 



-fe- % t) 2 /(2 CT |t) _ p 2 Vj z jd /^ e -(z j -2z jd -r 1j t) 2 /(2a]t) 



(20) 



The probability of company j surviving to time t is the integral of this over allowed values 
of zj and it is given by 



PAt) 



1 



erf 



2a]t 



+ 1 | - l e *ii*i*tf [ erf 



Vjt + z jd 



2a]t 



+ 1 



(21) 



dx e 3 ' 



where the error function is defined by 

2 f z 

erf ( z ) = ~m \ 
V 71 ^ Jo 

The leading joint survival probability is 

Pi 9) Jt) = P 1 (t)---P n (t). 



(22) 



(23) 



The contribution to the joint survival probability density that is linear in the asset 
correlations pW(zi, . . . , z n ) satisfies the differential equation 



d d 1 d 2 \ 1 

+ J2 - 2 S J P (1 \ Z ^ ■ ■ ■ i Z n, t) = - Wit 

i i / i=ij 



d 2 



dt 



13 dzidzj 



p(°\z 1 ,...,z n ,t) = 0. 

(24) 
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The Greens function (see, for example, [13] for a review of perturbation theory techniques) 
for the diffusion equation, K(z±, . . . , z n , xi, . . . , x n ; t — t'), satisfies the differential equation 



' d d 1 d 2 \ 

7u + H^/f - ~ oE^p K ( z ii ■ ■ -> z n,xi, • ..x n ,;t-t') = Sit-t^Sfa-xi) ■ ■■8{z n - 



8t ' ^' n dzi 2^" l 3zf 
i i 1 / 

The solution to the differential equation for p^\zi, . . . , z n , t) is 

^ roo 



- X 7i 



(25) 



roo roo Q2 

/ dxi---/ dx n K(z 1 ,...,z n ,x 1 ,...x n ,;t-t') - p (0) (xi, . . . , x n , t') (26) 

Jx ld Jx„ d OXiOXj 



The Greens function satisfying first passage boundary conditions is 

K(zi, . . . , z n ,x 1 , . . . x n , ;r) = 9{t)p{z 1 ,x 1 ;t) ■ ■ ■ p{z n , x n ; r), 
where 9(t) is defined to be unity for r > and zero otherwise and 



(27) 



p{Zj,Xj-,T) 



-(zj-Xj-Vjt) 2 /(2<t?t) _ ~(z :j +x :j -2z jd -r 1:j T) 2 /(2a 2 T)+2r 1j (z jd -x j )/. 



(28) 

Note that p(zj, Xj; r)dzj is the probability that the asset value of company j starts off at 
Xj and random walks without defaulting in the time r to a region dzj about the value Zj. 
Hence p(zj,0;r) = p(zj,r). Using Eqs. (27) and (28) 

p {1 \z 1 , . . . ,z n ,t) = j J[ Pk(z k ,t) \ J diloiOjiij 

roo Q \ / roo Q \ 

J dxip(zi,Xi;t - t')—p(xi,t')j [J dx j p{z j ,x j ;t-t')—p{xj,t')\ (29) 



We are interested in the effects of the asset correlations on the joint survival probability. 
Using the above results we find that 



^l„w/^l„(*) = 5E 



,(0) 



1 



1 



2^m)p 3 {t) jo 



dt'cjiCTjiijAij^t'), (30) 



where, 

Aij (t,t) 
and 



Xj d 



d 

dxiUi(xf,t- t')—pi(xi,t') 
dxi 



/ \ 

dxjUj{xj;t - t')—pi(xj,t') , (31) 
x jd ox j ) 



7]jT - Z jd + Zj 



2a]r 



+ 1 | - I e 2%(^ -%)/-| [ erf 



VjT + Z jd - Zj 



2a]r 



+ 1 



(32) 
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Eqs. (30), (31), and (32) are the main results of this paper. They hold even in the case of 
time dependent correlations with the replacement £jj — ► £ij(t') in Eq. (30). It is possible 

(2) 

to derive the expression for P± (i.e, the correction to the joint survival probability at 
second order in the correlations) however the number of numerical integrations required to 
compute it increases (from three) to six. 

Note that Eq. (30) implies that at linear order in the correlations the joint survival 
probability of n firms is expressible as a linear combination of the joint survival probability 
of pairs of firms. Explicitly, 



Pl,...,n( t )= P l( t )--- P n(t) 



1 



Pi At) 



(33) 



Thus, for example, Pi >2)3 (t) = A, 2 (i)P 3 (t) + Pi )3 (f)*M*) + ^.sftW*) " 2P 1 (t)P 2 (t)P 3 (t). 

It is sometimes convenient to introduce the random variables hi(t) that take the value 
1 if firm i survives without defaulting until time t and zero otherwise. Expectations of 
products of these variables are equal to the joint survival probabilities, E[n\(t) ■ ■ ■ h n (t)] = 
Pi,...,n{t)- Fluctuations about the expected value of these variables are characterized by, 
8hi(t) = hi{t) - E[hi(t)\. Eq. (33) implies that 



E[8hi(t)5hj(t)8h k (t)\ =0. i^j^k. 



(34) 



To understand the implications of Eq. (34) for portfolio risk consider a portfolio of corporate 
bonds subject to default risk and take the investment time horizon to be t. For simplicity 
we assume an initial investment of one dollar in the bonds of firm % has, after time t, the 
value 1 + Cj dollars if firm % doesn't default and value Ri + c, if firm i does default. Ri is the 
recovery fraction and a is the promised corporate return. Then the random variable that 
represents return for firm i bonds over the time horizon t is = Cj — (1 — Ri)(l — hi{t)). 
In this model the return fj only takes on two possible values and hence is extremely far 
from normally distributed. Other sources of risk which we have not included 4 smooth out 
the probability distribution for r$ . If the assets in this portfolio are uniformally distributed 
amongst n bonds then Eq. (34) implies that the portfolio return's skewness goes to zero as 
n — ► oo. To get a non-zero value for the portfolio return skewness in this limit one must 
include in the joint survival probabilities terms second order in the asset correlations 5 . 
The default correlation dik{t) between firms i and k is defined by 



dikit) 



Pi,k(t) - Pi(t)P k (t) 



~ Pi(t))Pi(t)y/(l ~ P k (t))P k (t) ' 



(35) 



At linear order in the asset correlations the joint survival probabilities can be expressed in 
terms of default correlations, 



Pl,...,n(t) = Pl(t) ■ ■ ■ Pn(t) 



l + i^^(t) V(1 



mm - m)) 



Pi(t)m) 



(36) 



4 For example, changes in firm i's credit quality which do not result in default. 

5 For a discussion of the importance of skewness for portfolio allocation to corporate bonds see [14]. 
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In the case where all the off diagonal elements of the correlation matrix are taken to be 
the same £y = £, for i / j we can define the joint survival probability correlation duration, 



Dl,...,n(t) 



_i dP!,..., n (t) 



(37) 



which characterizes the sensitivity of the joint survival probability to changes in the asset 
correlations. Our results give an explicit formula for the joint survival probability correlation 
duration evaluated at zero correlation, 



Consider five hypothetical firms with identical properties Oi = 0.30 (annualized), di/ Viq=0.30 
and £ij = £ , for i / j. Here we calculate their risk free joint survival probability and so 
we set Hi = r. Finally we take Aj = r so their default barriers grow at the risk free 
rate (we also assume no dividend payments, qi = 0). Then we find that Pj(5yr) = 87.3% 
and so pjty 3 4 5 (5yr) = Pj(5yr) 5 = 50.6%. The effects of correlations are calculated by 
performing the integrals in Eq. (30) and (31) numerically. The joint survival probabil- 
ity correlation duration is D± i 2,3,4,5 (5yr) = 10 x 0.30 2 x 0.611 = 0.55. This implies that 
5 (5jv)/P^ 5 (5yr) = 0.55 x £. So for an asset correlation £ = 0.30 the correlations 
increase the five year survival probability from 50.6% to 50.6 x (1 + 0.55 x 0.30)% = 58.9%. 
Default correlations increase with the time horizon t. For example, Lucas [15] estimates 
that over one year, two year and five year time horizons default correlations between Ba 
rated firms are 2%, 6% and 15% respectively. Using Eqs. (30), (31) and (32) we find that 
with £ = 0.3 the five year risk neutral default correlation for any pair of the hypothetical 
firms considered in this example is djj(5yr) = 11.3%. We will provide evidence in Section 
4 that £ = 0.3 is small enough that first order perturbation theory is an accurate approx- 
imation. Hence our results are useful for pricing first to default baskets of five firms with 
realistic values for the correlations. 

3 Asset Correlation Dependence of The Joint Survival Prob- 
ability For Five Industrials 

As a practical example of the results of the previous section we calculate here the depen- 
dence of the risk neutral joint survival probability of five basic industrials on their asset 
correlations. The five companies we consider are: Alcoa Inc. (ticker AA), Dow Chemical 
Company (ticker DOW), E.I. du Pont de Nemours and Company (ticker DD), International 
Paper Company (ticker IP) and Weyerhauser Company (ticker WY). 

For these companies we take di = £^,0/^0 where Dio is the present value of the total 
debt. This corresponds to a recovery fraction on the total debt very near unity. Note that 
the recovery fraction on the total debt can be much greater than on a particular bond issue 
since the recovery fraction on a particular bond issue depends on its level of subordination. 
Sometimes practitioners take the recovery fraction to be random [16] reflecting uncertainties 




■t 



Dl,..., n (t)k=0 



dt'(JicrjAij(t, t'). 



(38) 
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associated with accounting transparency and other issues . The initial asset value is V^o = 
Sio + Dio where Sio is the present market cap. Finally the corporate (annualized) asset 

value volatilities are determined from the historical stock volatility a\ (taken over a two 

(S) 

year period) via, Oi = (Sio/Vio) x a\ . Table I gives the input parameters for each of these 
companies that are needed to compute the dependence of their joint survival probability on 
their asset correlations. The last column of Table I gives their five year risk neutral default 
probabilities calculated using Eq. (21). Note that the dividend yield in the fourth column 
is with respect to the asset value and not the equity value. We calculate the risk neutral 
probabilities so the drifts \Xi are set equal to the risk free rate. Finally we assume the default 
thresholds grow at the risk free rate (i.e., \ = r) and hence the survival probabilities are 
independent of the value of the risk free rate. 

The small value for the default probability for E.I. du Pont de Nemours and Company 
is a consequence its very small default threshold. There may be items (e.g. legal liabilities 
) that effectively increase the value of di/V^ and are not reflected in our estimate. 

Table I: Input parameters and five year risk neutral default probabilities for five 

industrials. 



Ticker 


di/Vio 




Qi 


1 - ^(5yr) 


AA 


0.19 


31.2% 


1.5% 


4.7% 


DD 


0.089 


25.2% 


2.9% 


0.02% 


DOW 


0.24 


25.0% 


2.6% 


3.6% 


IP 


0.39 


16.5% 


1.4 % 


2.6% 


WY 


0.47 


16.5% 


1.4 % 


8.3% 



Our goal in this section is to calculate the dependence of the risk neutral five year joint 
survival probability on the asset correlations. Since all the companies are in the same sector 
we assume their asset correlations are the same £jj = £, for i ^ j. Then the joint survival 
probability correlation duration Di^.s^yr), evaluated at zero correlation, characterizes this 
dependence. Using the parameters in Table I and Eq.s (31), (32) and (38) we find that for 
the five companies considered here, 

£>i,...,5(5yr)|f=o = 0.036. (39) 

The five year joint survival probability in the uncorrelated case (i.e., £ = 0) is 0.82 and so 
the probability of one or more of the companies defaulting in a five year period is 

1 - A,..., 5 (5yr) = 0.18 - 0.029 x £. (40) 

Since we have worked to linear order in perturbation theory Eq. (40) receives corrections 
of order £ 2 . For £ = 0.30 correlations decrease the probability of one or more defaults 
happening in five years from 18% to 17%. The correlations have small effect because 
Weyerhauser Company has a much larger probability of defaulting in this period than the 
other firms. 

6 It is possible to generalize the results of this paper to the fluctuating default threshold case. Then 
Eq. (31) would contain two additional integrations over values of di and dj weighted with the probability 
distributions for these companies default thresholds. 
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4 Comparison With Multivariate Normal Copula Function 



In the practitioner literature a multivariate normal Copula function is sometimes used to 
determine the dependence of joint survival probabilities on asset correlations. A multivariate 
normal Copula functions gives for the joint survival probability of n companies 7 , 



1 /*oo roo 

A,...,n(t) = 77, W2 /XT? / dxi ■ ■ ■ / dx n exp 



~ (2vr)«/ 2 VdW J X1 (t) ^ " 7x4*) [ * 5 



(41) 



where, £>„■ , are elements of the inverse of the nxn asset correlation matrix. The thresholds 
are determined from the analogous formula for the survival probability of a single firm i. 

R t (t) = -= d Xi e~*i/ 2 (42) 

V27T J Xi (t) 

which gives 

Xi (t) = v / 2err 1 (l-2P i (t)) (43) 

where erf -1 denotes the inverse of the error function (i.e., erf _1 (erf (x)) = x) Expanding 
in the correlations we get for the case of a multivariate normal Copula function that the 
correlations cause a first order correction to the survival probability given by the simple 
expression, 

To study the difference between the value for survival probability with a multivariate 
normal Copula function a first passage model we consider a hypothetical situation where all 
the companies have the same properties (i.e., the volatilities, default thresholds and initial 
asset values, and hence the probabilities of default). Then we can write, in either the first 
passage or multivariate Copula function models 

= ^Y^(Mt) («) 

where as before we have set £jj = £ for i ^ j and A does not depend on the correlations. 



Table II: A(5yi) in a first passage model and using a multivariate normal Copula function 

for various choices of parameters. 





di/Vio 


P*(5yr) 


Xi(5yr) 


A {p (5jv)/af 


^c(5yr)/^ 2 


0.30 


0.20 


96.5% 


-1.81 


0.0697 


0.0717 


0.30 


0.30 


87.2% 


-1.14 


0.611 


0.636 


0.30 


0.40 


73.8% 


-0.636 


2.06 


2.17 


0.35 


0.20 


91.6% 


-1.38 


0.223 


0.231 


0.35 


0.30 


78.5% 


-0.789 


1.08 


1.13 


0.35 


0.40 


63.4% 


-0.343 


2.71 


2.87 



The Copula function contains more information than this. If the thresholds are evaluated at times ti, 
. . . t n then it gives the joint survival probability, when company i survives to time ti. Here we focus on the 
case where all the times are the same U = t. 
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We compute A for several examples using a multivariate normal Copula function and 
using a first passage model with correlated random walks and compare them. The results 
are presented in table II. The thresholds Xi m the multivariate normal Copula case are 
determined using Eq. (43) which sets the probability of survival for each company equal 
to what it is in the first passage model. In columns five and six of table II we present 
values of A (divided by the volatility squared) using the first passage model (Af p ) and the 
multivariate normal Copula function model (Ac)- In the first passage case we assume the 
default threshold grows at the risk free rate and that there are no dividend payments. Since 
we are computing the risk neutral joint probability of survival the drifts are set equal to 
the risk free rate 8 . The first two columns give the values of the annualized volatility and 
initial default threshold divided by asset value we use. Column three gives the resulting 
five year survival probability for a single company and column six gives the threshold in the 
multivariate normal Copula function that corresponds to this survival probability (computed 
using Eq. (43)). Note that the values of A agree to within about 10% in these two models 
and for the same correlations the multivariate normal Copula function gives a larger joint 
survival probability than the first passage model. The agreement between the multivariate 
normal Copula function and the first passage model is better for higher quality firms that 
have a larger survival probability. 

Since the multivariate Copula function allows us to compute the survival probability 
for any correlation we can use it to see over what range of asset correlations first order in 
perturbation theory applies. First consider the parameters used in the first line of table 
II (i.e., Xi = —1.8102). In that case we find that perturbation theory gives -Pi,...5(5yr) = 
0.837+0.054^. On the other hand explicit evaluation of Eq. (41) gives for £ = 0.1, 0.2, 0.3, 0.4 
and 0.5 the joint survival probabilities Pi,... 5 (5yr) = 0.842,0.849,0.858,0.867, and 0.877 
respectively. Since in this case -Pi,...5(5yr) is close to unity comparing linear perturbation 
theory with the value of 1 — -Pi,...5(5yr) is the best way to measure of its accuracy. For 
£ = 0.3 linear perturbation theory for the probability of one or more of the five companies 
defaulting during the five year time horizon is accurate to about 3% while for £ = 0.5 linear 
perturbation theory only accurate at the 8% level. Next consider the parameters used for 
the second line of table II (i.e., Xi = —1-1383). In that case we find that perturbation theory 
predicts Pi,...5(5yr) = 0.5056 + 0.289£ while explicit evaluation of Eq. (41) gives for £ = 
0.1, 0.2, 0.3, 0. 4 and 0.5 the joint survival probabilities Pi,... 5 (5yr) = 0.534, 0.563, 0.591, 0.619 
and 0.648 respectively. Evidently in this case first order perturbation theory is a reasonable 
approximation for the joint survival probability even if asset correlations are not small. 
For £ < 0.5 it is accurate to better than 1%. Linear perturbation theory for the joint 
survival probability works better for companies of lower credit quality that have higher 
default probabilities. 

5 Concluding Remarks 

We have studied the impact of correlations on joint survival probabilities using a structural 
first passage model. We derived, using perturbation theory in the correlations, an explicit 

8 Then the survival probabilities don't depend on the risk free rate and so we don't have to assume a 
value for it. 
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formula for the impact of asset correlations on the joint survival probability. We compared 
these results with what a multivariate normal Copula function predicts for joint survival 
probabilities and found (in the cases we considered) that the agreement between the two 
models was quite good for higher quality firms with a low default probability over the time 
horizon under consideration. Our results use first order perturbation theory which includes 
effects linear in the asset correlations. 

At first order in perturbation theory the joint survival probabilities can be expressed 
in terms of default correlations. It is possible to extend our analysis to second order in 
perturbation theory. 

We argued that for five firms first order perturbation theory provides a useful approxi- 
mation to the joint survival probability for realistic values of the asset correlations. When 
the number of firms is very large our results provide a limit where numerical simulations 
of correlated random walks, with first passage boundary conditions, can be compared with 
analytic results. 

One advantage of structural models is that companies bond and stock properties are 
related. Using historical stock volatilities we computed the dependence of the joint survival 
probability for the five basic industrials (i.e., Alcoa Inc., Dow Chemical Company, E.I. du 
Pont de Nemours and Company, International Paper Company and Weyerhauser Company) 
on their asset correlations. 

The methods we use are applicable when the asset correlations are time dependent. One 
possible extension of the work presented here is to a situation where asset correlations are 
typically small but for a short period of time they are large. As long as the time period 
where they are large is small compared with the total time horizon perturbation theory can 
be used. Periods where correlations are large could result, for example, from market stress 
associated with heightened geopolitical risk. 
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